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Abstract 

We focus on two dependency quantities of a max-stable random field X on some 
space T: the extremal coefficient function 9 which we define on finite sets of T and 
the extremal correlation function x( s ^) = hm^-foo ¥(X S > x \ X t > x). We fully 
characterize extremal coefficient functions 9 by a property called complete alterna- 
tion and construct a corresponding max-stable random field. Simple properties and 
consequences concerning the convex geometry of extremal coefficients are derived. 
We study how the continuity of X, 9 and x are linked to each other, and we show 
that extremal correlation functions x allow for convex combinations in general, and 
for products and pointwisc limits if the resulting function is continuous. These are 
operations which are well-known for positive definite functions, but the latter are 
non-trivial for extremal correlation functions. Finally, we regard some additional 
implications, when the random field X on T = R d is stationary. 
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1 Introduction 

Max-stable random fields have gained increasing attention in many areas of application 
such as hydrology or meteorology - see e.g. [J 0, [2H for some recent examples. The 
structure theory and classification of max-stable processes is well-studied, exemplarily we 
refer to [8, 14, 3Q]. However, appropriate modelling of the spatial dependence structure is 



still a difficult question. In contrast to Gaussian random fields, the dependence structure 
of a finite sample of a max-stable random field cannot be described parametricaliy. 
Instead, the notion of spectral measure, stable tail dependence function or dependency 
set (cf. (H, @, 0, [22I], for example) provide three equivalent concepts each of which 
fully describes the dependence structure of a multivariate max-stable distribution with 
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fixed marginals. Further, Pickands dependence function 2l[ comprises the complete 
dependence structure of a bivariate max-stable distribution. 

However, all of these dependency descriptors may have a very complex structure 
which makes them hard to be estimated from data (cf. (9I. Eol. 11] and references therein 
for some approaches on the issue in a multivariate setting). 

Therefore, several attempts have been made to propose meaningful, but simpler 
dependence measures: In this paper, we focus on extremal coefficients [29( in a spatial 
set ting and a function which we call extremal correlation function and was suggested 
in I26l | as an analogue to the correlation function for extreme values. The latter is also 



a special case of the so-called extremoqra m Ml and sometimes referred to as (upper) 
tail dependence coefficient (function) [3, 0, 0], x- measure 0; 0] or extremal coefficient 



function [13]. Besides being a bivariate dependence measure, the extremal correlation 
function also contains mixing information of the process H. An estimator which stems 
from a geostatistical approach is proposed in [20[j. 

The purpose of this paper is to characterize these spatially defined quantities - 
extremal coefficients and extremal correlation function - through their properties, to 
give necessary and sufficient criteria and to establish operations on the class of such 
quantities. 

The text is organized as follows: In section [2] we give the definitions and our notation 
of the basic objects that we study, i.e. max-stable random fields and related quantities 
such as extremal coefficient function and extremal correlation function are introduced. 
In section [3] we characterize extremal coefficient functions by a concept called complete 
alternation and construct a corresponding max-stable random field, thus generalizing 
results of j25| to a spatial setting. After deriving some immediate consequences from 
this, section H] is an intermezzo towards the convex geometry related to the results of 
section [3l In section [5] we analyse how the continuity of extremal coefficient function, 
extremal correlation function and corresponding max-stable random fields are mutually 
linked. Section [6] classifies the convex sets of extremal coefficient functions and extremal 
correlation functions within other related convex sets of functions. We shall see that 
some of these convex sets of positive definite functions are additionally closed under 
pointwise limits and products. Finally, we present some additional conclusions which 
apply to stationary max-stable random fields on IR d in section [7J 



2 Foundations and Definitions 

We consider max-stable random fiels X = {Xt}t£T on an arbitrary location space T with 
the same one-dimensional non-degenerate marginals. For convenience we use standard 
Frechet marg inals, i.e. F(X t < x) = e~ l/x for t G T and x > 0. Then, max-stability of 
the random field X simply means that for any n independent copies X«,...,X( n ) of 
X we have C{n~ l maxj = i i ... )n (XW)) ~ C{X), i.e. these fields are equally distributed. 

We will frequently fix a finite subset M C T and use the following conventions: 
Elements of M M , the vector space of real- valued functions on M, are denoted by x = 
(xt)teM where xt = x(t). By 1^ G IR M , we denote the indicator function of L where 
I / L C M. Thus, the elements {l¥,} teM G R M form a standard basis of R M . The 
standard scalar product on this space is (x,y) = J2teM x t ' Vt- Analogous notation 
is made for the space [0,oo] M , the space of [0, oo]-valued functions on M. In this 
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context, we define the expressions l/oo := and 1/0 := oo. Further, we consider some 
reference norm ||-|| on IR M (not necessarily the one from the scalar product) and denote 
S M ■= {a G [0,oo] A/ : ||a|| = 1}. 

The full dependence structure of each multivariate distribution of {Xt}t^M for some 
finite subset M C T can be described by using one of the following equivalent approaches: 



spectral measure Hm (cf. [a, 1221], for example), i.e. the measure on Sm such that 
for x G [0,oo] M 

-logPpf t <x tj *€M)= / max ( ^] H M {da) (1) 

Js M teM \ x tJ 



stable tail dependence function Im (cf. for example), i.e. the function on 

[0, oo] Af defined through 

£ M (x) := - log¥(X t < l/x t , t G M) = f max (a t ■ x t ) H M (da) 

J Sir * eM 



• dependency set Km (cf. [19l|). i.e. the unique compact convex set Km C [0, oo) A/ 
satisfying 

£m{x) = sup{(x,y) : y G Km}, (2) 
where it suffices to check this condition for ^ x G [0, oo) M . 

Of course, if we consider such quantities for different subsets M C T, they have 
to be consistent among each other to define a random field X = {Xt}t£T on T. For 
instance, due to our choice of standard Frechet marginals we have the normalization 
£ M (lf t} ) = I Sm a t H M (da) = 1 for t G M. 

For any non-empty finite set of locations M C T its extremal coefficient is defined 
using some x G (0, oo) by 

afi\*\ -logP(max teM ^i < x) f . . , , . M . , 

9{M) := _ lpgp(x ^ x) = max (a t ) H M (da) = i M (l M ). (3) 

Indeed, 6{M) does not depend on x G (0, oo) and takes values in the interval [1, \M\]. 
Roughly speaking, the extremal coefficient 9{M) grasps the effective number of inde- 
pendent variables among the random variables {Xt}teM- It is coherent to set 9{%) := 0. 
Consistency of the finite dimensional dependency measures requires also that 0{L) = 
£l(1V) = £m{1¥) for ^ L C M, for example. We will consider 9 in this paper as a 
function on the set of finite subsets ofT, which we denote by J-(T). The function 

9 : T(T) [0, oo) 

is then called extremal coefficient function. It is an application of PHopitals rule that 
9(M) coincides with 

a . 1 - P (maxtgM X t <x) P (3 1 G M such that X t > x) 
9{M) = hm — ; = hm — r . 4 

y ' xtoo 1 - P {X t < x) xtoo P (X t > x) K J 
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That is why the extremal correlation function 6 is intimately connected with a bivari- 
ate characteristic of the max-stable random field X which we call extremal correlation 
function \ given by the following limit 

X :T xT -> [0,1] x (s,t) : = lim P(X S > x \ X t > x) (5) 

xfoo 

This quantity dates back to [3] and has entered the literature under different names, 
most prominently (upper) tail dependence coefficient (function) [l|, Q, E3] or ^-measure 
[HIH]- We prefer the term extremal correlation function, sincey is a symmetric positive 
definite kernel on T x T as has been observed already in 0, 0, for example. The 
relation 

6({s,t}) + X (s,t) = 2 Vs,teT (6) 
holds because of @ and is widely known [l|, M, G3] • For technical reasons we introduce 
r, :TxT4 [0,1] r ] (s,t):=l- X (s,t) = 9({s,t})-l. (7) 
In terms of the spectral measure, we can express these bivariate functions as 

x (s,t) = min (a s , a t ) H M (da) r)(s,t) = - \a s - a t \ H M {da). 
J Sm •> Sm 

for any finite M CT with s,teM. 

3 Harmonic analysis of extremal coefficients 

3.1 A consistent max-linear model for max-stable random fields 

Suppose we are given a collection of real numbers (Tjf) where M C T ranges over all 
non-empty finite subsets of T and L over all non-empty subsets of M. We consider 
random fields X* = {X^}t^T on T with marginal distributions of the following form: 

- log F(XT <x t ,teM)= V t¥ max — VM € F(T) \ {%}, x G [0, ool M 

« *6i X t 

(8) 

Remember that J~(T) denotes the set of finite subsets of T. Such a model has been 
studied already in [251 ] for finite sets T. We draw our attention now to arbitrary index 
sets T and list the conditions on these numbers (r^f) which guarantee that they define 
the distribution of a max-stable random field X* with standard Frechet marginals: 

First of all observe that once these numbers (r^f) define a valid random field via 
model ([8]), max-stability is guaranteed immediately. As to the other conditions: 

• Formula ([8]) defines a valid distribution function for all non-empty finite M C T if 

t^>0 VM € F(T) \ {0}, i^LcM (9) 

In case these inequalities are satisfied, the distribution of {X£}t£M coincides with 
the distribution of the random vector {Xt}t£M which is given as a max-linear 
model built from 2' M I _1 i.i.d. standard Frechet variables {Vl}$^lcm by X t = 
maxteLcMT^VL (cf. 0). 
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• Subsequently, Kolmogorov consistency is ensured if and only if 



_ MU{t} MU{t} 
T L ~ T L + T LU{t} 

which is equivalent to 

T K — 2-*i T KUJ 
JCM\A 



VM G T{T) \ {0}, ®^LcM,teT,tg M, (10) 
VM G T{T) \ {0}, $ ^ K a A<z M. (11) 



The model has standard Frechet marginals if and only if 



r {t} 



1 



Vt G T. 



(12) 



To illustrate the model (|8j) better, we take a look at its spectral measure: Since it is 
a max- linear model, the spectral measure H* M of {X£ }t£M is discrete and can be written 
as 



17* II i M II c 

Q^LcM l|i£ J ll 

with respect to the reference norm ||-|| on M M (cf. {!])), see also figure [TJ 




12:i 
123 



Figure 1: Spectral measure representation of {Xl}teM for M = {1,2,3} if we 
choose the reference norm on M M to be the maximum norm. In this case the spec- 
tral measure simplifies to a sum of weighted pointmasses on the vertices of a cube: 
H 



M 



3.2 Complete alternation of extremal coefficients 



In [19|] and [25( the class of possible extremal coefficients has been fully characterized 
in a multivariate setting. This has been done by a set of inequalities which means the 
complete alternation of those coefficients and immediately transfers to the spatial setting 
as we shall see next. Behind the following characterizations, we use the fact that J~(T), 
the set of finite subsets of T, forms a semigroup with respect to the union operation U 
and with neutral element the empty set 0. The following notation is mainly adopted 
from and 0. For a function / : J-~(T) — > R and elements K,L G J-(T) we set: 



(A K f) (L) :=f(L)-f(LUK) 



5 



Definition 1. (complete alternation / negative definiteness) 

A function tp : F{T) — > M. is called completely alternating on J~(T) if for all n > 1, 
{K u ..., K n } C F{T) and K G T{T) 

(A Kl A K2 ...A Kn ^)(K)= Yl (-l)'^UuU^ <0. 

/C{l,...,n} V ie/ / 

A function ip : J-(T) — > R is called negative definite (in the semigroup sense) on T(T) if 
for all n > 2, {-ftTi, . . . , K n } C .F(T) and {a±, . . . , a n } C M with Yll=i a j = 1 

n n 

i=i fe=l 

Because the semigroup (T(T),U,$) is idempotent, these two terms coincide, i.e. i/j : 
J-{T) — > K is completely alternating if and only if ^ is negative definite (in the semigroup 
sense) (cf. 0, 4.1.16). 

An important example of a completely alternating function on J~{T) is the capacity 
functional C : J-(T) -> 1 of a binary random field Y. It is given by 



C(M) = P(3 f 6 M such that Y" t = 1) = P f max = 1^ . 
The following theorem characterizes such capacity functionals: 



(13) 



Theorem 2 ([la], p. 52, for example). The function C : T(T) — > E is £fte capacity of 
a binary {0, l}-field on T with the same one- dimensional marginals if and only if the 
following four conditions are satisfied: 

(i) C takes values in the interval [0,1]. 

(a) c(0) = o 

fraj C({t}) =: p < 1 is independent oft&T. 

(iv) C is completely alternating. 

Surprisingly, we may characterize the possible extremal coefficient functions on T(T) 
quite similarly: 

Theorem 3. a) The function 9 : J-(T) — > R is the extremal coefficient function of 
a max-stable random field on T with standard Frechet marginals if and only if the 
following three conditions are satisfied: 

(i) 9{%) = 

(ii) 0({t}) = 1 for allteT 
(Hi) 6 is completely alternating. 
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b) If these conditions are satisfied, the following choice of coefficients 
rf := -A {tl} . . . A {tl} 9(M \ L) = \ L) U I) 

VMeJ(T)\{0},fl^L = {ti,... ) t i }cM 

/or model defines a max-stable random field X* on T with standard Frechet 
marginals which realizes 6 as its own extremal coefficient function 0* . 

Proof. If is an extremal coefficient function of a max-stable random field X on T, then 
necessarily 0(0) = and 6({t}) = 1 for all i € T (cf. Q). Because of @ we have 

„, , „ , P (3 1 G M such that X t > x) , (M) 
0(M) = lim — ; r ~ ' = hm 

where denotes the capacity functional for the binary random field Yj = lx t >x and 
p( x ) = EYf = 1 — e _1//;c . Since complete alternation respects scaling and pointwise limits, 
complete alternation of follows from theorem[2j This shows the necessity of (i),(ii),(iii). 

Now, let : T(T) — > E be a function satisfying conditions (i),(ii),(iii) and let (r^f) 
as given above. We check that the numbers (rff) fulfill the (in) equalities (|9j) . (fTU|) , (fT2 |) . 
Indeed we have: 

© The inequalities = -A {tl} . . . A {tl} 6(M \ L) > follow directly from the 
complete alternation of 0. 

(gH) For t € T we have r { { * } } = {t} = 1. 

(|10p From the definition of Am we observe 

Cg } = - A m A {u • • • a w 0((m u {t» \ (l u {*})) 

= -A {tl} . . . A M 0(M \ L) + A {tl} . . . A {tl} 9{M U {t} \ L) 

_ -M Mu{t} 

Thus, (t^) define a max-stable random field X* on T with standard Frechet marginals 
as given by model ([8]). Finally, we compute the extremal coefficient function 0* of X* 
and see that it coincides with 0: 

**(M)§ £ rf= £ ^(-1)I^0((M\L)UI) 

Q^LcM ®^LcM ICL 

= E w E (-i) |xnL|+1 = E w(-(-i*=m)) = 0(m).d 

M\LCK 

So, theorem [3] is a complete characterization of all possible extremal coefficient func- 
tions : J-(T) — » [0, oo) of max-stable random fields, and it provides a max-stable 
random field X* corresponding to a given valid extremal coefficient function. In [2H] the 
last issue of the proof is derived for finite sets T by a Moebius inversion. The relation to 
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their proof becomes more transparent if we compute 9*(A) for A C M from {Jl^^lcM 
directly: 

nA) § £ £ £ r M j= £ r M (14) 

0^A- C A O^tfcA JCAI\A LcM-.LnA^tt 

Further, it follows that the bivariate distribution of the random field X* is given by 

-logP(X s * <x,X; <y) = r)(s,t) min [\, +max(^,-^) , (15) 
where rj = 1 — % (°f- ©)• We will use this later. 



x y J \x y 



3.3 Consequences of complete alternation of extremal coefficients 

Here, we collect some immediate consequences of theorem [3] that concern the closure 
of extremal coefficient functions 9 and extremal correlation functions x °f max-stable 
random fields, the operation of Bernstein functions on such functions, and a reinterpre- 
tation of the role of coefficients rff from model ([8]) in terms of an integral representation 
of 9. 

Closure Since 9 may be characterized by a list of (in)equalities as given by theorem 
O we may take pointwise limits (which may be read as "setwise" limits since the points 
of J-(T) are in fact sets): 

Corollary 4. The set of extremal coefficient functions 9 : J-{T) — > [0, oo) of max-stable 
random fields on a given space T is closed under pointwise convergence. 

Corollary 5. If T is countable, the set of extremal correlation functions \ '■ T x T — > 
[0, 1] of max-stable random fields on a given space T (cf.^)) is closed under pointwise 
convergence. 

Proof. Let Xn ■ T x T — > [1, 2] be extremal correlation functions of max-stable random 
fields for n 6 N, such that we have lim n _ i . 00 Xn(s, t) — > f(s,t) for all (s,t) G T x T. Let 
9 n : J-~(T) —7- [0, oo) be the corresponding extremal coefficient functions for n € N. On 
each finite set M, the value 9 n {M) is bounded by \M\. Thus, each 9 n in fact defines a 
point in the following product space V (by evaluation): 

9 n = {9 n {M)) M&HT) G {0} x J] [1, \M\] =: V 

MeJ 7 (T)\{0} 

The topology of pointwise convergence on .F(T) corresponds to the product topology on 
V C [0, oo)" F ( T ) . Since the intervals [1, \M\] are compact and the product is taken over the 
countable set J-(T) (since T is countable), the sequence of extremal coefficient functions 
9 n : F(T) — > [0, co) possesses a convergent subsequence 9 n i — > 9 with respect to pointwise 
convergence. Corollary 0] now tells us that 9 is again the extremal coefficient function 
of a max-stable random field on T. The observation that f(s,t) = lim n f^. 00 Xn'( s T't) = 
lim n /^. 0O 2 — 9({s, t}) = 2 — 9({s, t}) shows that the limit function / on T x T is again 
an extremal correlation function of a max-stable random field. □ 
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Operation of Bernstein functions Since condition (iii) in theorem[3]is equivalent to 
6 being negative definite, we get that Bernstein functions operate on extremal coefficient 
functions. 

Definition 6. (Bernstein function) 

A function g : [0, oo) — >■ [0, oo) is called a Bernstein function if one of the following 
equivalent conditions is satisfied (cf. 4.4.3 and p. 141) 



(i) The function g is of the following form for some c, b > and v a positive Radon 

lo i+a' 



measure on (0, oo), such that J* °° j^u(d\) < oo: 



g (r) = c + br + J (l- e _Ar ) u(d\) 



(ii) g is continuous and g € C°°((0, oo)) with g > and (-l)^^ 1 ) > for all n > 0. 
(Here, g^ denotes the n-th derivative of g.) 

(iii) g is continuous, g > and g is negative definite as a function on the semigroup 
([0,oo),+,0). 

For a comprehensive treatise on Bernstein functions including a table of examples 
see [23[ ] . Since Bernstein functions operate on negative definite kernels (cf. [3], 3.2.9 and 
4.4.3), we have the following corollary: 

Corollary 7. Let 9 : J-(T) — > R be an extremal coefficient function of a max-stable 
random field and g a Bernstein function which is not constant. Then the function 

m „ mm - <y> 

9(1) - 9(0) 

is again an extremal coefficient function of a max-stable random field. 

For instance, if 9 is a valid extremal coefficient function of a max-stable random 
field , then also log(l + 0)/log(2) or 9 q for < q < 1 are valid. Of course, this leads 
also to the corresponding transformations on the set of bivariate quantities x an d f] (cf. 
©,©)• 

Further, the characterization of theorem [3] allows to derive many simple inequalities 
that an extremal coefficient function 9 of a max-stable random field necessarily satisfies: 

Corollary 8. Let 9 : J- — > R be an extremal coefficient function of a max-stable random 
field and g a Bernstein function. Then the function g{9 — 1) : J 7 \ {0} — > R satisfies the 
triangle inequality 

g (9{A UB)-l)<g (9(A uC)-l)+g (9{C UB)-1) for A,B,Ce ?(T) \ {0}. 

Proof. Since 9 is an extremal coefficient function of a max-stable random field, it is 
completely alternating which is equivalent to being negative definite (in the semigroup 
sense, cf. theorem [3]). Subtracting 1 does not change this property. But notice further 
that 9 takes values in {0} U [l,oo), where the value is only attained for the empty 
set 0. Thus, 9 — 1 : T \ {0} — s- R is negative definite and takes values only in [0, 00). 
Applying a Bernstein function g does not change this property. By 0] 8.2.7, this also 
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means that / := g{9 — 1) : J- \ {0} — > M. is completely alternating on T \ {0}. Since we 
have also / > on T \ {0}, we may derive for A, B,C € T \ {0} 

f(AuB)-f(AuC)-f(CUB) 
< /(C) + /(^ U B) - f(A U C) - f(C U B) 

= (/(C) - /(A U C) - f(C UB) + f(A U B U C)) + {f {A U B) - f(A U B U C)) 
= A A A B /(C) + A c /(iUB)<0. □ 

By letting yl = {s}, B = {t}, C = {r} in corollary [8j we may also derive such 
inequalities for the respective bivariate quantities: 

Corollary 9. The function r] = 1 — x °f a niax-stable random field on T (cf. ^) 
satisfies for any Bernstein function g the triangle inequality 

g o n(s, t) < g o n(s, r) + g o n(r, t) for s,t,r G T. 

Integral representation Theorem 2 also allows to reinterpret the role of the coeffi- 
cients ijf of the max-linear model ([8]) within the framework of integral representations 
of negative definite functions. 

Definition 10. A bounded semicharacter p on J-(T) is a function p : F{T) -4- {0, 1} 
that satisfies p(0) = 1 and p(A U B) = p(A)p{B) for A,B £ T(T). 

The set of bounded semicharacters on J-"(T) is denoted J-{T) and forms itself a 
semigroup with respect to pointwise multiplication with neutral element the constant 
function 1. 

Since we have theorem [3l we may apply 0] 4.1.16, that is any extremal coefficient 
function 9 of a max-stable random field allows for an integral representation: 

Corollary 11. Let 9 : J-(T) — > IR be an extremal coefficient function of a max-stable 
random field. Then 9 uniquely determines a positive Radon measure p on J-"(T) \ {1} 
such that 

9(A) = p({pef(T)\{l}\p(A) = 0}) 

where 9({t}) = 1 fort E T. The function 9 is bounded if and only if p(J r (T)\{l}) < oo. 

If we restrict ourselves to a finite non-empty subset McT, the restricted extremal 
coefficient function 9m ■ J~(M) — > [0, |M|] is still completely alternating and corollary [TT1 
applies, of course. But now, since we consider only a finite M, we can use the additional 
isomorphism of finite (topologically discrete) idempotent abelian semigroups: 

(JXM),-,1)-(^(M),U,0) p^ p| (M\L) 

Lep-Hi) 

in order to reformulate corollary [11] for this situation: Any such restricted extremal 
coefficient function 6m ■ F{M) — > IR uniquely determines a positive Radon measure pu 
on T(M) \ {0} such that 

e M {A) = p M {{LeT\{%} \AnL^®})= Yl ^(W) 

where X^teL MAf ({-^}) = 1 101 t £T. A comparison with (]14|) reveals that 

»M({L}) = rf. 
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4 Convex geometry of extremal coefficients 

We will fix again a finite non-empty subset M C T. Remember that we denoted by Km 
the dependency set of {Xf}teM ( CI *- ©) which was introduced in \\\$\ and which is the 
unique compact convex set Km C [0,oo) m satisfying 



hl(x) = sup{(x,y) : y G K M } VO + x G [0,oo) M . 

Here, we will show that the dependency sets K* M of our max-linear model (jHj) play an 
exceptional role among all dependency sets sharing the same extremal coefficients: 

Theorem 12. Consider a max-stable random field with extremal coefficient function 9 
and denote Km its dependency set for a finite non-empty M C T . Further denote 

H L :={xG [0,oc) A/ : (x, if ) < 9{L)} 

the half space in [0, oo) bounded by the plane with defining equation (x, if ) = 0(L). 

a) Then K M C Do^lcJI/^- 

b) For each ^ L C M there exists a point x L in the intersection x L G Km H 

c) The dependency set K* M of (X*)t£M from model |2p with extremal coefficient function 
9 is given by K* M = {\^ LcM Hl- 

Proof, a) Let ^ L C M and x G K M - Assuming (x, if) > 9{L) we would have the 
contradiction 9{L) = £ M (lf ) = sup{(x, if) : x G K M \ > 9{L). So (x, if) < 9{L) 
and part[aj) is proven. 

b) Since Km is compact, we have 9(L) = ^(lf) = sup{(x, if) : x G Km} = 
(x(L), if) for some x L G Km- This shows part lb]). 

c) Denote the right-hand side by 

C M ■= P| H L = {xe[0, oo) M : (x, if) < 9(L) for all ^ L C M}. 

0^LcM 

The inclusion /C^ C £m is proven in part[aj). So, we still need to show the other 
inclusion Cm C /Cj^. By definition (cf. ([2])), /C^ C [0, 1] A/ is the unique compact 
convex set satisfying £m{x) = sup{(x, y) : y G i^^fl for 7^ x G [0, oo) M . The closed 
convex may also be described as the following intersection of half spaces (cf. 



271 ] . section 1.7.): 



Km= fl [0,oo) M : (x,y)<£ M (x)}. 

xes M 

Thus, it suffices to show the following implication in order to prove Cm C K* m : 
x G S M and y G £ M => (x, y) < Ijif(i) 
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We now prove this implication: Without loss of generality, we may label the elements 
of M = {t 1; . . . , t m } such that x tl > x t2 > ■ ■ ■ > x tm . Then we may write x € Sm C 
[0,oo) M as 



X = + fon-l ~ 1 M\{t m } +■■■ + ( X t2 ~ gfaj , l{ tl ,fa} + (ftl - ^fe) jfa} 

>o >o >o >0 

Taking the scalar product with y E £tv/ ; we conclude 
<x J i/><x tm e(M) + (x tn _ 1 -x tm )e(M\{t m }) + ... 

• • • + (x i2 - x t3 )e({t!, t 2 }) + (x tl - x t2 )e({tx}) 

= x tm (8(M) - 6(M \ {t m })) + ■■■+ x t M{h, t 2 }) - edh})) + xtMh}) 

(16) 

On the other hand the stable tail dependence function £m of model (|8j) is by this 
ordering of the components of x given as 



e M (x)= t l 



M 



max St 
teL 

Q^LcM i=l \LcM:ti,...,U-i£L,tieL 



From formula (|14|) we see that this expression of £m(x) coincides with the right-hand 
side of (fT6|) . Thus, we have our desired inequality (x,y) < £m{x). This shows part 
Ej). □ 
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Figure 2: Dependency set /C M of {X?} teM of model © for M = {1, 2, 3}. 

So, if we fix the extremal coefficient function 6 of a max-stable random field on T, 
then model (jH|) yields a maximal dependency set K* M w.r.t. inclusion, that is 



K 



M 



U 



/Cm- 



/Cm dependency set 
with the same extremal coefficients as /Ct 
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It is an open problem and it would be interesting to know whether there exist also 
minimal dependency sets and if they would help to better understand the classification 
of all dependency structures. A very naive idea would be to take for each L in each of 
the sets {x E [0, oo) M : (x, 1%) = 9(A)} n Dlcm^^l for ^ A C M one point and 
then to take the convex hull with included. However this fails to be a dependency set 
in dimensions \M\ > 3, since it is not even a zonoid, which would be necessary (cf. jl9l|). 

5 Continuity 

Continuity of functions on J-(T) In this section we require T to be a metric space. 
We technically prepare the notion of continuity that we will use in connection with 
extremal coefficient functions 9 : J-(T) — > [0, oo). Therefore let / : J-(T) — > R be a 
function on the finite subsets of T. Then / induces a familiy of functions {/^}m>o 
where f( m ' : T m — > M. is given by 

f^\t 1 ,...,t m ) = f({t 1 ,...,t m }). (17) 

We see that for each m the induced function /( m ) : T m — > R is symmetric, that is for 
any permutation tt 6 S m 

..,t m ) = /M(t 7r(1) ,...,t 7r(m) ). (18) 

Further, the family /( m ) : T m — > M. is consistent in the following way: For any non- 
negative integers ki, ■ ■ ■ , k r with Yli=i h = m and (ti, . . . , t r ) € T r we have 

/^(t^^tl, . . . ,tr 1 _^tr) = f( r \h, ...,t r ). (19) 

ki times k r times 

On the other hand, given functions {f^ m '}m>o which satisfy (|18p and (|19p . they define 
a well-defined function / : J-(T) — > R via (|17p . 

Definition 13. Let / : J-(T) H-Mbea function on the finite subsets of a metric space 
T. We say that / is continuous if all induced functions : T m — >• R are continuous 
for all m > 0. 

Remark 14. To be more precise, this notion of continuity for functions on J-(T) used 
in this definition corresponds to the final topology on J-(T) when we view it as the 
direct limit J-(T) = lim J-" m (T) where J rm (T) denotes the set of finite subsets of T with 
number of elements less than or equal to m, and where J^l^T) = T m / ~ carries the 
quotient topology of T m . 

Continuity of 9 and \ For an extremal coefficient function 9 : J-(T) — > R we denote 
its induced functions as described in the paragraph above by 0( m '. 

Theorem 15. Let X be a max-stable random field on a metric space T with standard 
Frechet marginals, 9 its extremal coefficient function and x ^ s extremal correlation 
function. Further, let X* be the random field from model |5|) that realizes 9 as its 
extremal coefficient function, and thus also x as ^ s extremal correlation function (cf. 
theorem^). Then the following implications hold: 
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a) X is stochastically continuous. ==? 9 is continuous. ==? \ ^ s continuous. 

b) The function 77 = 1 — x ( c f- satisfies for any e > 

F(\X* S - X*\ > s) < 2 - exp (-^)) < f »?(*,*)• 

x is continuous. X* is stochastically continuous. 

Proof, a) The second implication is obvious from x( s i~k) = 2 — 9^ 2 \s,t). As to the first 
implication: Stochastic continuity of X means that for any e > 0, for any t € T 
and sequence -> t we have P(|X t o) — Xt| > e) — > 0. From this, we can easily 
derive that for any e > and any (ti, . . . , t m ) G T m and a sequence (i^ \ ■ ■ ■ ; im^) ^ 
(ti,... , i m ), also P(||(Z(„) — 2^)2=1 1| > e) — >■ for any reference norm ||-|| on M m . 
The latter implies the corresponding convergence in distribution: F ' („) <„) — > 

l^t-^ ,...,6772 J 

F(t!,...,t m )- Since logF( tlj tm ) : [0, oo) m — > E is monotone and homogeneous, we 
have that for x > the point (x, . . . , x) € (0, oo) m is a continuity point of F(ti,...,t m ) 
(cf. [22I], p. 277). Thus, we may conclude that the induced function 6^ on T m is 
continuous, since 8^ m \ti, ...,t m ) = -x\ogF( tli ... jtm )(x, ... ,x). 

b) Let e > 0. We will prove the statement for 2e instead of e. Therefore, consider 
the following disjoint events on a corresponding probability space (0, A, P) for k = 
0,1,2,... 

A k := {lo € n : X t *(u/» G (fee, (fc + 2)ef \ ((k + l)e, (k + 2)e] 2 } 

The disjoint union {J^LqA^ is a subset of {uj £ Q, : \X*(u) — X£(lu)\ < 2e} and so 

(00 \ 00 n 

\jA k )= j2 F ( A k) = j™Z) p (^*)- 
fc=0 / k=0 n °° k=0 

For further calculations we abbreviate for p, q € N U {0} 

B(p,q) :=F(X* S <p-e,X; <q-e). 

Note from (fT5|) that B(p,q) = B(q,p) and £?(p, 0) = 0. With this notation we 
rearrange 

n n 

P(^fc) = -B(n + 1, n + 1) + 2^ [B(k + 2, k + 1) - B(k + 2, fe)] 

fc=0 fc=0 
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For the second summand we have (cf. (|15|) ) 

n 

[B{k + 2, k + 1) - B{k + 2, k)] 



k=0 

k=0 

n 

= ^exp 

k=0 
n 

> ^exp 

fc=0 



exp 



k + 2 

T](s,t 



+ 



1 



k + 2 

y(s,t 
2e 



k + l 
exp 



exp 



r](s,t) 1 
k + 2 k 



1 



(k + l)e 



exp 



1 

ke 



exp 



1 



(fc + l)e 



exp 



2e 



exp 



1 



(n + l)e 



Finally, 



X;| > 2e) = 1 - P(|X S * - X t *| < 2e) < 1 - lim Vp(4) 



fc=0 



1 + lim B(n + 1, n + 1) - 2 lim V [J3(ife + 2, jfe + 1) - B(k + 2, k)] 



< 1 + lim exp 

n— >oo 



" (n + l)e 



fc=0 
2 lim ( exp 

n— >oc 



2e 



exp 



1 



(n + l)e 



2 exp 



2e ~ 2e' y 1 1 



c) This follows from part |b|) . 

The following corollaries are immediate consequences of theorem [T5l 



□ 



Corollary 16. Let : J-(T) — > [0, oo) be the extremal coefficient function and x '■ 
T x T — > [0, 1] t/ie extremal correlation function of a max-stable random field on a 
metric space T with standard Frechet marginals. Then 



is continuous. 
Corollary 17. Let T be a metric space. 



X is continuous. 



a) The set of continuous extremal coefficient functions 9 : F{T) — > [0, oo) of max- 
stable random fields on T coincides with the set of extremal coefficient functions 
6 : J-(T) — > [0, oo) of stochastically continuous max-stable random fields on T. 

b) The set of continuous extremal correlation functions x '■ T x T — > [0,1] of max- 
stable random fields on T coincides with the set of extremal correlation functions 
X ■ T x T —?■ [0, 1] of stochastically continuous max-stable random fields on T. 

Finally, we will show that we are allowed to take pointwise limits on the set of 
extremal correlation functions of max-stable random fields as long as the limit function 
is continuous: 
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Theorem 18. Let T be a separable metric space and let Xn '■ T x T —> [0, 1] be the 

extremal correlation functions of some max-stable random fields on T such that the 
functions Xn converge pointwise to a continuous function f : T x T — > [0,1]. Then f is 
an extremal correlation function of a stochastically continuous, max-stable random field 
on T. 

Remark 19. However, we are not aware of whether it is still true, that the limit function 
/ is an extremal correlation function of a max-stable random field on T if we drop the 
requirement of / being continuous. 

Proof of theorem [7#1 Let Q C T be a countable dense set; then Q m C T m is also 
countable and dense for any m > 0. Since Xn converge pointwise to /, the restrictions 
Xti\q x q also converge pointwise to /|q x q- Because of corollary [5] and the arguments 
in the corresponding proof, the restricted function f\n x n is an extremal correlation 
function \ of a max-stable random field X on Q, which is a random field of type ([8]). 
Let us denote its extremal coefficient function by 9. 

The idea of this proof is now to show that the induced functions 9^ : Q m — > [1, m] 
can be uniquely extended to continuous functions on T m . For m = and m = 1 this 
is trivial, since 6>(°) = and 0« = 1 are constant by definition. Secondly, since / is 
continuous on T x T the function 2 — f : T x T — > [1,2] is the unique continuous extension 
of9W = 2-f\ QxQ toTxT. 

Moreover, set r/ = 1 — /. Then r] is also continuous with n(t, t) = for all t € T. Since 
rj(s,t) is the pointwise limit of functions which are symmetric and satisfy the triangle 
inequality (cf . corollary [9]) , the function r\ itself is symmetric and satisfies the triangle 
inequality r](s,t) < r](s,r) + i](r,t) for all s,t,r G T. 

Now, fix m > 3 and consider 9^ : Q m — > [1, m]. Note that even though we can 
derive continuity of 9^ m ' from corollary [To] we may not derive a continuous extension 
directly. But the following property of 9^ will be enough: Let B$(t) C T m denote the 
open (5-ball around t 6 T m . 

Property A For t = (t\, . . . , t m ) € T m and e > we can find 5 > such that for any 
q,qe B s (t)nQ m we have \9^(q) - 9^ m \q)\ < e. 

This means, that #( m ) : Q m C T m — > [1, m] is Cauchy continuous. In fact, because of 
the symmetry, it suffices to establish the following property: 

Property B For t = (t±,... ,t m ) € T m and e > we can find 5 > such that for any 
q,qe B s (t) n Q m we have 9^ m \q) - 9^{q) < e. 

Interchanging the roles of q and q in property B yields property A with the smaller 5 
from property B. Now, we prove property B in three steps: 

1st step For all E\ > there is <5i > 0, such that for all q, q € B$ 1 (t) PI Q m we have 
2m T,iLiV(Qi,Qi) < El- 
Proof of 1st step: Since s i— > i](ti,s) is continuous at ti 6 T for i = l,...,m, there 
exist op > such that r](ti,qi) < for q, L € BM)(ti). Choose 5\ > such that 
Bsj^it) C YYh=i B m (U) . Let g,g £ S^^). Then the triangle inequality of r\ yields: 

mm m 

2m viQuQi) <2m^2 + viU, Qi)] < 2m ^ 2^—^ = e 1 

i—l i=l i=l 
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2nd step For q G Q m denote F q (x) := F(X qi < x,...,X qm < x). For all e x > 
there is Si > 0, such that for all q,q € Bs 1 (t) fl Q m and for all a,x > we have 
F q (x) <F q (x + a) + ^. 

Proof of 2nd step: A priori we have for any a, x > and all q,q S Q" 1 

m 

F g (x) = F(X qi <x,...,X qm <x,Y, \X qi -X qi \ < a) 

8=1 

m 

+ P(X gi <i,...,4 <x,J2\X qi -X 4i \ >a) 

8=1 

= P(X ?1 < x,. . . ,X ?m < x, \X qi -X qi \<a,..., \X Qm -X q J< a) 
+ F(\X qi -X^\>^or ... or \X qm -X q J>^ 

<F q (x + a) + J2^(\ x qi - x ^\ > ^) <^K* + a ) + ^&(*>«) 

8=1 8=1 

The last inequality follows from theorem [T5l part lb|) . Now, the 1st step gives a suitable 
choice for 5i for given E\. 

3rd step Proof of property B: 

We may assume e < 2. Choose a > such that 2am < e. In particular this implies 
e~ m > e~^ and we may choose e\ such that on the one hand s\ < a{e~ m — e~2°) and 
on the other hand E\ < ae _m /(| — 1). This ensures: 

(- log (e- -fl))<f a»d p^^if (20) 

a 

Finally, choose 5 := 5i as in the 2nd step (or equivalently 1st step) for this £\. Then we 
have for arbitrary q, q € B$(t) that 



M (q) ~ {m) {q) = (1 + a) (- log F q (l + a)) + log F q (l) 

< (1 + a) (- log (X(l) - ^) ) + log F,(l) 



a 

logF,(l) - log (F q (l) - ^)) + a (-log - ^ 

Here, the first equality is just the definition of 8^ and the second inequality due to the 
choice of 5 (cf. 2nd step). Now, we apply the mean value theorem to achieve 

§W(q)-6( m \q) = max \ ■ — + a (- log ( FJl) - — 
Ce[F,(i)-a.,F fl (i)] £ a V V a 

.fi + o(-log(F a (l)-^ 



Using e~ m < F 9 (l) and ([20]) . we arrive at 
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This proves property B. 

Thus, for each m > we may uniquely extend 9^ to a continuous function #( m ) : 
rpm _^ [i )772 ,]. Naturally, this extension {#( m )} m >o inherits the properties (JTSJ and (JT^D 
from {9^} m >o- Thus, it defines a function on F(T) via (|17p . By construction, this 
extension respects the three properties of theorem [3l i.e. 9 defines indeed the extremal 
coefficient function of a max-stable random field X* on T as constructed in theorem 
[3l As mentioned above, / = % is the corresponding extremal correlation function x 
of this random field. Because of theorem [T5l the random field X* is stochastically 
continuous. □ 



6 Inclusions of convex sets of functions 

6.1 Extremal coefficient functions 

We consider the following sets of functions on F(T) in this section: 

logP(max ie M X t <x) X a max-stable r.f. on T 



6 max (T) := M4 9{M) 



log F(X t < x) 



with std. Frechet marginals 



e Um (T) :-- 



M H> lim 



logP(maxtgju Z t < x) 
logP(Z t < x) 



Z a r.f. on T with the same 
one-dimensional marginals 
and upper endpoint oo, 
such that the limits exist 



e bin (T) := 



M i — y 



(max tgM Y t = 1) 
F(Y t = 1) 



Y a r.f. on T with the same 
one-dimensional marginals 
taking values in {0, 1} ; 
such that Elf + 



In particular, we denote m ax(^) the set of extremal coefficient functions of max- 
stable random fields on T, which is already characterized by theorem [3] as the set of 
completely alternating functions 9 on F{T) with 9{%) = and 6({t}) = 1 for t 6 T. 

Proposition 20. The following inclusions hold: 

a) closure(9 bin (r)) = Q max (T) = 6 lim (T) 

b) e bin (T) = {9 G G max (T) : 9 bounded } 

where the closure is meant with respect to pointwise convergence. 

Proof, a) Firstly, let F 6 0bin(^). Then F is a normalized capacity functional (cf. 
(|13p ). such that F{%) = and F({t}) = 1. Therefore, F is completely alternating 
(cf. theorem [2]) and satisfies all conditions of theorem [3l i.e. we have F S ©max(7")- 
This shows 



e bin (T) c e max (T) 



(21) 
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Secondly, let F G Oii m (T) for some random field Z. Similarly to flD) it is an easy 
exercise to show that 



F(M) - 1' 1 ~ F ( maXfgA/ Zt - x ^> - r IP (max f sM t Zt >x = 1) 
[ l-F{Z t <x) F{t Zt > x = l) 

Therefore, F G closure(0bi n (7 1 ))- This shows 

ei im (T) c closure(G bin (T)) (22) 

Additionally, we have by definition that ©max (7") C 9iim (T) and ©max CO is closed 
(cf. corollary U|) . So, all inclusions follow from 

|22} |2U 

©max(T) C 0ii m (T) C closure(©bi n (r)) C closure(G max (T)) = 6 max (T). 

b) Let F G ®bin{T) and Y a corresponding random field. Then F is bounded by Let 
G 6 max (T) be bounded, say by K. Clearly, K > 6({t}) = 1. Set C(A) := 0(A)/K. 
Then C satisfies all requirements of theorem [2] and defines a binary field Y as in the 
definition of @ hin (T), such that P(max 4eM Y t = l)/F(Y t = 1) = 0(M). □ 

In fact, all involved sets are convex: 

Proposition 21. The sets ©bin(^) an d 8 max (T) = 0ii m (T) are convex. 

Proof, a) ©bin(r) is convex: Let Yp and Yq be binary random fields on Tcorresponding 
to F,G G 9 bin (T) and p F := E(Y F ) t > 0, p G := E(Y G ) t > 0. Choose £/ ~ Bin(l,a), 
independent of Yp and Yq. Then t/Yp + (1 — U)Yq is binary and its capacity 
C a : .F(T) — > K is given in terms of the capacities Cjr and by 

QPF ^ (1 - a)p G 

Gq = — T, s — l f H — r: v — <-,g- 

apF + (1 - aJ^G apF + (1 - aJpG 

Now, for every (3 G (0, 1) there is a unique a G (0, 1) such that ap ^L a ) p = ft- 

b) ©max(r) is convex: Consider the the max-combination X = max(aJi, (1 — a)^) of 
two independent random fields X\ and X2 on T with standard Frechet marginals and 
extremal coefficient functions 9\ and 02, respectively. Then the extremal coefficient 
function 9 of X is a9\ + (1 — a)02- □ 

Remark 22. Alternatively, the convexity ©max (7") = 0ii m (T) follows directly from theo- 
rem or we could argue that it must be convex as the closure of the convex set ©bin (J 1 ). 
But the argument given here is a bit more constructive. 
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6.2 Extremal correlation functions 



Similarly to the previous section, we define the following sets of functions on T x T: 

X a max-stable r.f. on T 



H max (T) := <^ (s,t) ' y X (s,t) = lim F(X S > x \ X t > x) : 

xtoo with std. Frechet marginals 



Z a r.f. on T with same 

.a i- m,/^ ^ i ^ ^ \ one-dimensional marginals 
■Aim(T) := < (s,t) i-> lim P {Z s > x \ Z t > x : . 6 

xtoo with upper endpomt oo, 

such that the limits exist 



(s,t)^F(Y s = l\ Y t = l) 



Y a r.f. on T with same 
one-dimensional marginals 
taking values in {0, 1}, 
such that Elf + 



All of these are sets consist of symmetric, positive definite functions onTxT with 
values in [0, 1], which take the value 1 on the diagonal. 

Proposition 23. The following inclusions hold: 

a) S bin (T) C H max (T) C H lim (r) 

b) Si im (T) = closure(H bin (T)) = closure(H max (T)) 

where the closure is meant with respect to pointwise convergence. 
If T is a separable metric space, we have additionally 

c) {x € H max (T) : x continuous} = {x€ Hn m (T) : x continuous}. 

Remark 24. In particular if T is countable, we have E max (T) = Eii m (T) (cf. corollary [5]) . 
It would be interesting to know whether also in the general situation H max (T) = Hi; m (T) 
is correct. 

Proof of proposition \23l a) Let {Ytjt^r be a binary random field with same one-dimensional 
marginals and EY t ^ 0. Then ¥(Y S = 1 | Y t = 1) = 2 - C ({s , t}) / C ({t}) in terms 
of the capacity functional C of Y, just as x( s >t) = 2 — 8({s,t}) for the quantities 
that we defined for max-stable random fields (cf. ([6])). Therefore, the assignment 
C I—?- 9 from capacity functionals of binary random fields to extremal correlation func- 
tions of max-stable random fields, which is given by 6{M) := C(M)/C({t}) proves 
Sbin(r) C E max (T) (cf. the proof of Qhm(T) C 6 max (T) in proposition [20]) . The 
inclusion E max (T) C Hi; m (T) follows directly from their definitions. 

b) Because of[aJ it suffices to show Hii m (T) = closure(Hbi n (? 1 )). By definition of Hn m (T) 
and Sbi n (T) and considering the random fields Y t = lz t >x indexed by x > 0, we have 
^iim(T') C closure(Hbin(7")). To show the reverse inclusion let f n —¥ f be a pointwise 
converging sequence with f n G E^ in (T) and Y n the corresponding binary random 
fields, which are chosen to be independent. For further computations we abbreviate 
a n {s,t) := F({Y n ) s = l,(Y n ) t = 1) and b n := P((Y n ) t = 1) = E(Y„) f (which is 
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independent from t). Note that by choice of Y n we have < a n (s,t) < b n < 1 and 
< b n and that ""^'^ — >. £) a s n — > oo. We may assume that one of the following 
two cases applies: 

First case: There exists a lower bound b > such that b n > b for all n € N. 
Second case: The sequence b n converves monotonously to as n — y oo. 

(Otherwise there exists a subsequence (b nk ) — > that converges monotonously to as 
k — > oo and we work with the sequence (fn k )k instead of (f n )n which still converges 
to / pointwise.) 

Set X = N ■ Yjv, where N is an independent random variable taking values in N with 
p n ■= p(JV = n) = 2~ 2 " for n > 2 and p x := F(N = 1) = 1 - £ 4 = 2 2 " 2 '- Then the 
sequence (p n ) nS M satisfies the property that Y^= m +i — ~~ * as m — > oo. 
Now, 

P(X S > x | X t > x) = 11 



for m= \x~\. In order to prove lim^oo F(X S > x \ X t > x) = f(s,t), it suffices to 
show that lim m _>. 00 Y"^™ , i ^- • — = 0, since < a n (s,t) < b n . Therefore, let us 
consider the two cases: 

First case: Then £~ m + i fe " £ < \ En= m+ i XT= m +i •£ by 

choice of (p n )neN- 

Second case: Then £~ m+1 fc • £ < E~= m+l 1 " ^ ^ by choice of (p n ) ne ®. 
c) This is a direct consequence of theorem 1181 □ 

Again, all involved sets are convex: 
Proposition 25. All sets S max (T), Eu n (T), Si; m (T) ore convex. 

Proof. The convexity of Hbin(r) and H max (T) follow from proposition EH since we have 
F(Y S = 1 | Y t = 1) = 2 - P(max(Y s , F) = l)/P(Tt = 1) for the binary random field 
F and x(a,t) = 2 - 6»({s,t}) (cf. ©). The convexity of H lim (T) = cl(E hin (T)) = 
closure (H max (T)) follows from the convexity of E^i n (T) or E max (T), since the closure 
does not affect the property of being convex. □ 

Further operations on the sets H max (T), Efr m (T), En m (T) 

Corollary 26. Let f E Hii m (T). Then 1 — / satisfies for any Bernstein function g the 
triangle inequality g(l - f(s, t)) < g(l - f(s, r)) + g(l - f(r, t)). 

Proof. This is a consequence of proposition [23] and corollary [9j □ 

Corollary 27. a) Let fi and be functions both either from E^^T) orE\[ m (T). Then 
the product f\ ■ fi belongs again to this set of functions. 
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b) Let T be separable metric. Let fi and fi be functions from 

{x G H max (T) : x continuous} = {xPO G S max (T) : X stochastically continuous}. 
Then the product f\ ■ $2 belongs again to this set of functions. 

Proof, a) The statement is easy to prove for Hbm(T') by taking the product of random 
fields of independent random fields. For Hji m (T) it then follows from proposition 1231 

b) Since E max (T) C closure(Hbi n (T')) (cf. proposition [23]) we may choose sequences 
/n,lj/n,2 G Hbm(r) which converge pointwise to f\ and /2, respectively. Now, for 
each n the product /„ 5 i • f n ^ is again an element of Ebin(r) C H max (T) and pointwise 
/n,i • /n,2 converges to /i • /2- Now, the claim follows from theorem [HI Q 

Remark 28. In particular the set Hn m (T) is closed under convex combinations, products 
and pointwise limits, whereas the set 

{x G H max (T) : x continuous} = {x{X) G H max (T) : X stochastically continuous} 

is closed under convex combinations, products and pointwise limits as long as the limit 
function is continuous (cf. theorem I18f) . This allows for further operations built from 
these. We may apply convergent power series with positive coefficients: For instance, 
if x ^ Smax(T') is continuous, then also 1 — (1 — x) q G ^max(r) for q E (0,1). Note 
that corollary [7] is in accordance with these operations, since corollary [7J implies that 
for any x € H max (T) and any Bernstein function g with g(0) = and g(l) = 1, the 
function 2 — g(2 — x) 

is again in ^ max (T 1 ). However, the above is slightly more general 
for continuous members of H max (T). For instance, we get that 1 — g(l — x) € H max (T), if 
X £ "max(r) is continuous. Finally we want to remark that, although these operations 
(convex combinations, products, pointwise limits) are well-known operations on the set 
of positive definite functions, it has not been obvious so far that they are compatible 
with the restriction to (continuous) extremal correlation functions of max-stable random 
fields. 

7 The stationary case 

Most of the results stated before, are compatible with group actions on the location 
space T and we could formulate a version which incorporates this group action. For 
example if we let T = M. d and focus on stationary random fields on theorem [3] has 
the following version: 

Theorem 29. a) The function 6 : J-(T) — > R is the extremal coefficient function of 
a stationary max-stable random field on M d with standard Frechet marginals if and 
only if the following four conditions are satisfied: 

(i) 0(0) = 

(%%) 9{{t}) = 1 for allt G R d 
(Hi) 6 is completely alternating. 

(iv) 6(M + t) = 9(M) for any M G T(R d ) \ {0} and any t G R d . 
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b) If these conditions are satisfied, the following choice of coefficients 
rf := -A {tl} . . . A {tl} 9(M \ L) = ^(-1)^+^((M \ L) U I) 

ICL 

VM G.F(T)\{0}, 0/L = {h,..., ti} CM 

for model $8$ defines a stationary max-stable random field X* on R d with standard 
Frechet marginals which realizes 9 as its own extremal coefficient function 9*. 

Similarly, the results for the bivariate quantities transfer, in particular theorem [181 
and corollary [271 which leads to the respective operations on the set of extremal corre- 
lation functions of stationary max-stable random fields. 

Further, the triangle inequality 

g o r/(s ± i) < g o r](s) + g o r/(t) for s,t £ R d 

of the function rj = 1 — x of a stationary max-stable random field on R d and any 
Bernstein function g (cf. corollary [9|) can be found already in the literature for some 
specific functions g: The three choices g{x) = log(l + x), g(x) = (1 + x) T — 1 for 
< r < 1 and g(x) = 1 — (1 + x) T for r < as Bernstein functions yield precisely the 
set of inequalities derived in 0] Proposition 4. Therein a madogram approach is used 
instead. 

This triangular inequality where g = id is also the reason, why such functions 77 (and 
equivalently x = 1 — v) m &y not be differentiable at the origin unless they are constant 
[2a|. This is also true if we considered instead of R d any other abelian connected Lie 
group acting by translation on itself. For more implications of the triangular inequality 
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Finally, we want to mention that if we consider only stationary random fields on 
some of the inclusions in Proposition [20] and [23|) will be proper: 

Proposition 30. The following inclusions are proper inclusions: 



X (Y) G E hin (R d ) : Y stationaryj C | X (X) G H max (M d ) : X stationaryj (23) 

9{Y) G H b in(M d ) : Y stationary} C {(9(X) G E max (R d ) : X stationary} 

Proof. It suffices to show that the inclusion (|23p is proper. We will show that the l.h.s. 
of (|23p may not contain functions with compact support, whereas the r.h.s. does contain 
such functions: 

Suppose the l.h.s. contains a function / : R d ->■ [0,1], f(t) = P(Y t = 1 | Y = 1) 
with compact support for some stationary random field Y on ~R d with values in {0, 1} 
and p := EYo > 0. Then the normalized covariance function C(h) := Cov(Yj, Yo)/.P = 
/(f) - p has to be positive definite on R d . Further, set -y(t) := C(0) - C{t) = 1 - f(t), 
which is constantly 1 away from some compact set around the origin because / has 



compact support. Due to [15(], theorem 6, we conclude that 1 > 1 — p = C(0) > 
lim s _ i . 00 (2s)~ d Jj_ s s j d 7(£)df = 1, which leads to a contradiction. 

On the other hand, the r.h.s. of f)23|) contains functions of compact support. For 
instance, let A C R d be a compact subset of R d and let p = f l J 4(s)ds. Let II be a 
Poisson point process on (0, 00) x R d with intensity p~ 1 u~ 2 duds. Then the random field 
X t = sup( u >s ) e n ulA(t — s) on R d is stationary and max-stable [i^] and has extremal 
correlation function x(t) = p^ 1 f 1a{s — t)tA(s)ds, which has compact support. □ 
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